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Abstract Applying Gibb's geometrical methods to the thermodynamics of H-plasmas we explore the land- 
scape of the free energy as a function of the degrees of ionization and dissociation. Several approximations 
for the free energy are discussed. We show that in the region of partial ionization/dissociation the quantum 
Debye-Hiickel approximation (QDHA) yields a rather good but still simple representation which allows to 
include magnetic field and fluctuation effects. By using relations of Onsager-Landau-type the probability 
of fluctuations and ionization/dissociation processes are described. We show that the degrees of ioniza- 
tion/dissociation are probabilistic quantities which are subject to a relatively large dispersion. Magnetic 
field effects are studied. 
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1 Introduction 



In plasmas and in particular in hydrogen systems, ion- 
ization and dissociation processes play an important role. 
We find the equilibrium values of the degrees of ioniza- 
tion and dissociation by minimization of the free energy 
at fixed total mass P^[21S] • Further the ionization and dis- 
sociation rates are derived from kinetic equations 
Here we will start with a study of the geometry of the 
free energy landscape lElB] • The idea to apply geometrical 
methods in thermodynamical research goes back to pio- 
neering papers of Gibbs [7]. We apply Gibb's method here 
to plasmas. We will study different approximations to the 
thermodynamic functions and the influence of several ef- 
fects, our approach is limited to the region T > lOOOOi^T. 
In a preceding work we have developed a similar approach 
for ionization processes described by the Saha equation [6] . 
We include here also effects of magnetic fields which are of 
increasing interest not only for astrophysical applications 
but also for plasmas confined in magnetic traps. A prelim- 
inary work on magnetic field effects was published earlier 
[8]. The magnetic field introduces an anisotropy into the 
system. This changes the ideal gas contribution, the bound 
state energies and the contributions due to interactions of 
the free particles. Here we will take into account the ideal 
gas corrections in an exact way, and the bound state and 
interaction corrections in approximations which include 
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the quadratic orders in the field 0(i?^). Practically this 
restricts our approach to fields B < lO^T. 

We start with exploring the geometry of the free energy 
landscape as a function of the degrees of ionization and 
dissociation. On this basis we are able to describe not 
only the chemical equilibria but also fiuctuations around 
the equilibrium values [5]. 



2 The free energy density of hydrogen 

2.1 Defining the basic variables of the system 

We consider in the following hydrogen plasmas at fixed 
temperature T, and proton density n. Our aim is to study 
ionization processes 



H <=;> p + e^ 
and dissociation processes 

H2^H + H. 



(1) 



(2) 



For simplicity, the formation of H2 and species will 
be neglected. 

In order to define the thermodynamics of the system 
we introduce the free energy 



F (T, V, {N}) , 



(3) 



the density of the free energy 
F 



f {T,n) = — = f {T,ne,n,,na,nrn) , (4) 
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the degrees of ionization and dissociation [3] . We will use 
here the following definitions: 



a 



Pa = 



ni+na + 2n,; 



/3 



(5) 



We note that Pa is the relative amount of protons bound 
in atoms and Pm the relative amount bound in molecules. 
Due to the balance relation for the total proton density 

n, + na + 2nm- (6) 

we find the relations 

Pa = il-a)p, Pra = {l-a){l~P). (7) 

The condition of neutrality requires that electron and ion 
densities are always equal: 



(8) 



Therefore the free energy density depends only on 4 inde- 
pendent parameters 



/^/(T,n,a,/?), 



(9) 



If we include magnetic field effects we have 5 independent 
parameters 

f^f{T,n,B,a,p). (10) 

The equilibrium composition is found by the minimization 
procedure 

Instead of P we may use Pa or Pm as variational parame- 
ters. The relation 



a + /3a + /3m = 1 



(12) 



gives the possibility for visualization of the free energy 
surface on a simplex. 



2.2 The different contributions to the free energy 

We assume the following structure [2l[3|: 

f = + + fa + fl^ + fee + fu + he+ hs- (13) 

The first 3 terms describe the ideal contributions of the 
electrons, ions, atoms and molecules. The next 3 terms 
represent the Coulomb interactions electron-electron, ion- 
ion, and ion-electron. For these terms different approxi- 
mations exists as e.g. the quantum Debye-Hiickcl-approx- 
imation (QDHA) and Pade approximation p]. Basically 
we will use here the QDHA: 



fee ^" fa ^" fie ~ /. 



DH- 



(14) 



In the QDHA screening and quantum effects are both 
taken into account in a first approximation. The inclusion 
of further effects as Hartree-Focketc. by means of Fade ap- 
proximations makes no principal difficulties but increases 
the numerical efforts. We will show that in the region of 
partial ionization and dissociation the QDHA yields a rea- 
sonable description of ionization/dissociation phenomena. 
The interactions of the neutral components are taken into 
account by hard core approximation including effects of 
reduced volume [US]. 

We will discuss now these contributions in more detail: 

— the ideal free energy of electrons We use Fermi- 
Dirac statistics including excluded volume corrections 
representing the free energy in the following way [2] : 



fe 



, ^ ,'1 n,Al 



2 1 



(15) 



The function z(x) which is due to Zimmcrmann [2] in- 
terpolates between the low density limit: z{x) ~ Inx— 1 
and the high density limit: z{x) (x X3 . Further A^. = 
h I \f2jvm^kBT is the thermal wave length of the elec- 
trons, correspondingly we will use next the thermal 
wave lengths of ions, atoms and molecules (/li, Aa-, 
A^n). The excluded volume factor 770 < 1 expresses 
that some part of the total volume V which is occu- 
pied by atoms and molecules is not accessible to the 
electrons. The accessible volume (1 — 770)^ is smaller 
than the total volume and the effective densities are 
higher. A strict theory for ryo is not yet available. We 
will use here 



^7ri?^(na + 2n,; 



(16) 



as the definition for the packing parameter with an 
estimate for the atomic radius i?o = 0.78A. Effects of 
a magnetic field are included by the replacement [8] 



.A 



3tanh(a;e) 



(17) 



where Xe = hu!^/2kBT and uj^ = eB/m,,. 

the ideal free energy of the ions f^'^ which is given by 

the classical term 

= n.ksT [In {n,Al) - l] - ?i.fcBTln(l - 7]o). (18) 

The second term describes the excluded volume effect 
for the classical bare protons [2|. The parameter Xi = 
hujl/2kBT {ujI = eB/rrii) is very small Xi 1 for 
small B. Therefore the magnetic field effects can be 
neglected for protons. 

the ideal atom contribution: /^'*. We assume Boltz- 
mann statistics i.e. 



Pa'' = nakeT 



In 



^aAl 



- 1 



(19) 



including the internal states by a Brillouin-Planck-Lar- 
kin partition function [2] 



1 



1 



(20) 
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where ^ = 2\/Ry/fcsT with the ionization energy of 
hydrogen Ry ~ 13.598eV. Including the magnetic field 
we have to replace the atomic partition function CTq by 
an eflFective partition function, 



(^aiT) CTeft(T,B), 



(21) 



which reads [8] 



aeff(r,i?) = a,(ee) + g 



(22) 

with the interaction parameter = e^/47reo^_BT'Ae 
(Ae = h/ y/lmJiBT is the thermal wave length for the 
relative motion of electrons and protons (with an in- 
finit proton mass)) and [8] 



n 



2 



1 



1 



2! V2 

ei^Yl^ - 1 



n=l 



1 



3! I 2 



1 



2! \2j n4 
(23) 



the molecular contribution: We will use the repre- 
sentation [3] 



[inf"™^") 


- 1 







+/r+/r (24) 



including the internal molecular and atomic states 

'4.746eV\ 



(Tm = (Ja CXp 



(25) 



and the contributions of vibrational and rotational sta- 
tes [9] 



f 



vib 



= 6210K 



In 1 - e T 



n„,kfiT 



In- 



J. 7- J. 7- 



2T 
1 /T, 



T 3r 90 V T 



(26) 



(27) 



= 85K. 



the Debye-Hiickel contribution foH- Using here the 
so-called -^-approximation p] we arrive at the formula 



foH = -kBT——T {k a(T)) . 

iZTT 

Here the inverse Debye-radius is defined by 



' 2nie2 
eoksT 



(28) 



(29) 



and the effective thermal electron radius by 



(30) 



where Ag is the electron thermal wave length defined 
above. Further the r-function has the standard Debye- 
Hiickel form pj : 



t{x) = — 



ln(l 



X 

~2 



(31) 



The advantage of the QDHA is that magnetic field 
effects may easily taken into account by a factor [8] 



a(T) a{T,Xe) = a(T) 



1 



1 

48' 



(32) 



The magnetic field reduces the effective diameter. This 
reflects the fact that the magnetic field localizes the 
particles perpendicular to the field, 
the hard-core contribution fhs- For this term we use 
the Carnahan-Starling approximation [2] : 



fhs 



{ug + 2n^) (4?? - 3??^) 
(1-^)' 



(33) 



The volume fraction of the neutrals is here defined by 
the effective packing parameter 



47r 

V = y -R ("a + 2n„0 . 



(34) 



Following the work of Juranek and Rcdmer |10] we fix 
the effective packing radius of atoms at a mean value 
R = 0.37A. This is of course a rough approximation 
which might be refined by taking into account different 
radii for the atoms and for the molecules and by in- 
troducing temperature-dependent radii based on fluid 
variational theory pIJ] . 



3 Discussion of the approximations and 
influence of parameters 

For this section we neglect the magnetic field, i.e. we as- 
sume B = 0. We are interested in partial ionization and 
dissociation, i.e. in the region of the density-temperature 
plane where electrons, ions, atoms, and molecules exist. 
In this region the electrons may show quantum effects but 
the ions, atoms, and molecules are still classical. 

The excluded volume effect, i.e. the effect that the elec- 
trons and ions cannot penetrate into the neutral particles, 
are taken into account in the ideal contributions. The ex- 
cluded volume factor rjQ is defined above with a fixed ra- 
dius Rq. If this effect is not considered we find a lower 
ionization rate a at higher densities, shown in figure H 

A small difference of the radius i?o has an influence 
on the degree of ionization and dissociation, see flgure 
For larger radii, a begins to increase at smaller densities 



4 



Werner Ebeling et al.: Thermodynamics of hydrogen plasmas 



p[gcm 3] 
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10^^ 10^'^ 10''° 10""^ lO"""" 10'^ 10""- 10""^ 10""^ 
n [cm""^] 



Figure 1. Degree of Ionization and dissociation for different 
radii Ro at a temperature of T = 2 • lO'^K. From right to left: 
Ro = 0, Ro = 0.73A, Ro = 0.78A, and Ro = 0.83A 



Q [g cm ^] 

10~^ 10~^ 10~^ 10~^ 10~^ 10° 10^ 10- 




j_q18 -j^qW ■j_q20 ■j_q21 ^g22 ^q23 ■j_q24 ^q25 -j^q. 



Figure 2. Degree of Ionization as a function of the total 
proton density for different temperatures: (1) T = lO'^K, (2) 
T = 2 • 10*K, (3) T = 3 • 10*K, (4) T = 5 • 10'*K, (5) T = lO^K. 
The solid line represents the Calculation with the QDHA and 
the dashed line with Fade approximation. 

since the effective volume VcS = (1 — ilo)^, which is avail- 
able for the electrons, decreases with increasing tjq and 
thus increases the effective electron density nf^ = rie/Vcs. 
Therefore the system tends to higher degrees of ionization. 
The transition between low and complete ionization takes 
place in a very narrow density range. For high densities 
the excluded volume effect has to be considered and has 
still to be refined maybe by taking into account temper- 
ature dependencies of the radius [10] and by introducing 
occupation probabilities pJ^HHS]- 

The nonideal effects of the free charged particles can 
be described more correctly by means of Fade approxima- 
tion (see e. g. p3lLL5lQS|)- Figure [2] shows the ionization 
and dissociation rate as a function of the proton density 
calculated by using the QDHA and by the Fade formu- 
lae given in pj]. For the density range of 10^^ — lO^^cm"^ 




Figure 3. Free energy density at T = 300K, n = 2 • lO^^cm"^ 
represented on the simplex a, /3a, (3m- The minimum of the 
free energy is at the corner a — 0, f3a = 0, Pm = 1, i.e. in 
equilibrium we find practically only molecules. 

and the temperature range 1 • 10^ — 5 • lO^K we find only 
small deviations of the QDHA from the Fade approach. 
For higher temperatures above lO^K the deviations are 
much larger and the QDHA is less accurate. In the follow- 
ing we will restrict our study to temperatures much below 
lO^K, therefore it is justified to use further the QDHA to 
simplify the formulae and to reduce calculation efforts. 



4 The geometry of the free energy landscape 

On the basis of the formulae given in section [2] we repre- 
sented the geometry with respect to / (a, (3) at fixed T, 
n and starting with i? = (figures [SJ S}. We used the 
variables a, Pa, and j3m which allow a representation on 
a simplex according to eq. (p2J • The minimum of the free 
energy corresponds to the stationary degrees of ionization 
and dissociation. The form of the landscape around the 
minimum determines the thermal fluctuations. 

Depending on the values of temperature and proton 
density wc find different locations of the minimum. At 
T = 300K, ri = 2 • lO^^cm"^ wc observe a minimum in 
the corner a = 0, /3a = 0, (3m = 1, meaning that only 
molecules exist (figure S)- At different temperature and 
proton density, e.g. T = 4 • lO'^K, n = lO^^cm"^ we ob- 
serve a minimum in the center of the simplex. This cor- 
responds to a situation where free electrons and ions and 
bound states (atoms and molecules) exist under the same 
conditions (figure 3]). 

With the procedure described above we calculated the 
stationary degrees of ionization and dissociation for a lat- 
tice of points in the tcmperature-dcnsity-planc. In this way 
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Figure 4. Free energy density at T = 4 ■ 10*K, n = lO^^cm"^. 
The minimum is found in a region where electrons, ions, atoms 
and molecules are present. 




Figure 5. The degree of ionization a together with the disper- 
sion of Of as a function of the proton density and the tempera- 
ture. The coding is such that light regions correspond to small 
dispersion and dark regions correspond to large dispersion. 



we obtained the degree of ionization as a function of den- 
sity and temperature (figure [5) . 

We observe full ionization (a = 1 on a plateau in the 
upper right corner and at the left edge of the diagram) 
at very high densities and in the region of lower densities 
but high temperatures. At sufficiently low density we find 
only bound states when the temperature is low (a = in 
a region beginning at the lower edge of the diagram) . In 
the intermediate region of proton densities we find a val- 
ley of bound states stretching over the whole temperature 
range with small ionization taking place only at very high 
temperatures p]. 

Now we will look at the formation of molecules (figure 
[6) . We observe similar features as above in the correspond- 
ing parameter regions. At very high densities (left edge) 




Figure 6. The fraction of protons bound in molecules jSrr 
(1- a) (1-/3) 



no molecules exist as well as in the region of lower den- 
sities and high temperatures (upper right corner). When 
the density is sufficiently low and the temperature is low 
as well we find that the bound states observed in figure \5\ 
are practically only made up of molecules (/?„ w 1 at the 
lower edge). In the intermediate density region we find 
a mountain of molecules (/3m = 1) and a decline of the 
degree of dissociation with rising temperature. 



5 Fluctuations according to Boltzmann-Gibbs 

Let us first consider the general schema of the theory of 
isothermal fluctuations at given temperature T f[7W8\ . 
We follow Klimontovich who developed the general fluc- 
tuation theory of Boltzmann and Gibbs pS|- In the gen- 
eral case we study a set of intrinsic parameters A(X, a) = 
{Ai{X, a), . . . , An{X, a)) depending on dynamical variables 
X and external parameters a and define a conditional free 
energy F{a,T\A) which is the free energy at fixed val- 
ues A = A{X, a). The intrinsic parameters A{x, a) do not 
need to have any thermodynamic functions which corre- 
spond to them. The equilibrium value of the free energy 
is F(a, T). Then according to the Boltzmann-Gibbs prin- 
ciple the probabilities of fluctuations of A are given by 



p{A\a, T) ~ exp 



F{a,T) - F{a,T\A) 
kuT 



(35) 



In our ease the dynamical variables are A = {a, /3}, for F 
we choose the free energy per proton (j) = f/n with the 
equilibrium value (j)Q = fg/n at fixed temperature T. 

Then according to the Boltzmann-Gibbs principle the 
probability distribution of isothermal fluctuations reads 



p{a,f3\n,T) = 



(36) 



An illustration of the probability distribution calculated 
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a" 




e [gem 



Figure 7. The probabihty distribution in the a, /3-plane at 
T = 3 • lO^'K and n = 10^^ cm^ 



for a lattice of 50x50 points on the a, /3— plane is shown 
in figure \7\ 

Further we show in figure [5] the degree of ionization 
together with the dispersion of the degree of ionization 
coded in gray levels. The diagram was obtained with the 
following procedure: For a given temperature and proton 
density the stationary degrees of ionization and dissocia- 
tion (a, /3) were calculated. Then we additionally fixed the 
degree of dissociation at its stationary value [3 = [3 and 
calculated the probability p{a) to find arbitrary values for 
the degree of ionization. The dispersion of the degree of 
ionization {5a)'^ ~ (a — a)^ was then calculated as the 
root of the mean quadratic deviation from the stationary 
value a. 

We observe low dispersion (light) at very high densi- 
ties due to the hard-sphere contribution. High dispersion 
(dark) is observed in two transition regions, one when ion- 
ization processes start taking place at intermediate den- 
sities and temperatures (just above a = 0), the other at 
the transition to complete ionization (just below a = 1) 
where the remaining bound states break up. 

Inside the simplex we assume that the probability dis- 
tribution function p{a, /? = /?) is an one dimensional gaus- 
sian distribution for values a inside the borders and for 
small fluctuations Sa = a — a 



p{a) = C exp 



1 



2kBT V9a2 / 



(37) 



((5a) 2 = kBT 



(38) 



where a is the equilibrium value and therefore the most 
probable value. The constant C is determined from the 
normalization. The dispersion in the fluctuations is given 

by 

We see that the second deviation of 4>{a,f3 = (3), which 
is the curvation of the free energy landscape determines 
the dispersion {Sa)'^. Near to the corner of the simplex the 
distribution is highly non-gaussian and highly unsymmct- 
rical. In this case the border of frequent fluctuations (the 
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10- 



10- 



10" 



10^ 



10^ 



0.8 



0.6 



0.4 



0.2 




2q18 j_q19 j_q20 j_q21 j_q22 j_q23 j_q24 j_q25 j_q26 

n [cm~^] 

Figure 8. Degree of ionization as a function of the proton 
density n with the dispersion of the fluctuations of a at T = 
4 ■ 10**K. 



range of the dispersion) may be estimated by finding the 
roots of p{a) ~ p{a = a)/ ^/e. For gaussian distributions 
this criterion agrees with the one given above. 

In density-temperature regions where the minimum of 
the free energy is developed very well e.g. if the minimum 
is located in a corner of the simplex (see e.g. figure [3]) we 
have small dispersions. Is the free energy landscape flat 
e.g. if the minimum is located inside the simplex (see e.g. 
figure 3]) we have large dispersions. 

To further illustrate this, we show the dispersion of the 
degree of ionization at fixed temperature, corresponding 
to a cut through the diagram in figure [5] parallel to the n- 
axis. The dispersion is indicated by error bars, which are 
cut at the physically impossible values a > 1 and a < 
(figure [8}. Again the dispersion is small at high and low 
densities and is high near the transition a = 1 =^ a < 1 
where bound states emerge. 



6 Magnetic field effects 

We want to discuss now the influence of a constant uniform 
magnetic field on the ionization and dissociation equilib- 
rium. We have to carry out the replacements pT]), pT), 
and (321). These formulae equations are valid for weak 
magnetic fields, i.e. 



48 



< 1. 



(39) 



In figure Wi the degree of ionization and dissociation for 
various magnetic field strengths are plotted over the tem- 
perature. We find at a fixed density that the degree of 
ionization of a plasma in a magnetic field is higher com- 
pared to the field free case and increases with the field 
strength. For temperatures higher than lO^K and weak 
magnetic fields there are no differences in the ionization 
degree. Further we see in figure Qj]] that the curves devi- 
ate in the mid-density range where atoms and molecules 
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lO" 2-10" 3-10^ S-IO" 7-10" 10^ 

r[K] 

Figure 9. Degree of ionization and dissociation for different 
magnetic field strengths at a proton density n = lO'^'^cm"^. 



g [gem ^] 

io~^ 10'" io~^ icr^ io~^ 10° 10^ 10- 




10^^ 10^^ 10^^ 10^^ 10^" 10^^ 10^^ 



n [cm ^] 

Figure 10. Degree of Ionization and dissociation for different 
magnetic fields at a temperature T = 3 ■ lO'^K. 

exist up to the density where the transition to complete 
ionization takes place. 

7 Conclusions 

Based on a detailed study of the free energy landscape we 
developed a new description of ionization and dissociation 
phenomena. The relatively simple approximation QDHA 
used for the interaction term in the free energy allows the 
study of the dispersion of the degrees of ionization and dis- 
sociation and of magnetic field effects. The QDHA may be 
replaced by more refined approximations as Pade approxi- 
mation. We have shown that the study of free energy land- 
scapes in combination with the fluctuation theory yields 
relevant information on dispersion of the degrees of ion- 
ization and dissociation. The natural dispersion is usually 
in the range of 5 — 10% except on the borders. This has 
consequences for many experimental quantities as e.g. con- 
ductivity and optical properties. Since the conductivity is 
in first approximation proportional to a, i.e. a ^ aero, we 



may conclude that the thermal dispersion of the plasma 
conductivity a has nearly the same shape as the disper- 
sion of a. We note that the thermal fluctuations of a, (3, 
(7 etc. are in general not of Gaussian character what is 
due to the existence of borders as e.g. < a < 1. Study- 
ing the influence of magnetic fields we have shown that 
a remarkable influence on the degree of ionization starts 
at i? ~ lO^T, the effects increases with the magnetic field 
strength. At fields B > 5 ■ lO^T the degree of ionization 
and also the conductivity and other related quantities are 
substantially higher than for the field free case, the present 
approach however is no more valid at these extremely high 
fields. 

Recent experimental work on dense hydrogen concentrates 
mainly on Hugoniots p2] . The approach presented here is 
restricted to temperatures T > lOOOOif, therefore it can- 
not be used for the calculation of the low-temperature 
parts of Hugoniots. An extension taking into account re- 
cent results on the low temperature EOS of dense hydro- 
gen [2] is in progress. 

The authors thank Stefan Hilbert, Hauke Juranek, Ronald 
Redmer, Gerd Ropke, Manfred Schlanges and Werner Stolz- 
mann for many helpful discussions. The authors W. Ebeling 
and H. Hache profited very much from a stay at the Rostock 
University. 
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